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Motivated by recent experiments that reveal expansive fractional quantum Hall states in the n = 1 
graphene Landau level and suggest a nontrivial role of the spin degree of freedom [Amet et ah, Nat. 
Commun. 6, 5838 (2014)], we perform accurate quantitative study of the the competition between 
fractional quantum Hall states with different spin polarizations in the n — 1 graphene Landau level. 

We find that the fractional quantum Hall effect is well described in terms of composite fermions, 
but the spin physics is qualitatively different from that in the n = 0 Landau level. In particular, for 
the states at filling factors v — sf (2s =b 1), s positive integer, a combination of exact diagonalization 
and the composite fermion theory shows that the ground state is fully spin polarized and supports a 
robust spin wave mode even in the limit of vanishing Zeeman coupling. Thus, even though composite 
fermions are formed, a mean field description that treats them as weakly interacting particles breaks 
down, and the exchange interaction between them is strong enough to cause a qualitative change in 
the behavior by inducing full spin polarization. We also verify that the fully spin polarized composite 
fermion Fermi sea has lower energy than the paired Pfaffian state at the relevant half fillings in the 
n = 1 graphene Landau level, indicating an absence of composite-fermion pairing at half filling in 
the n = 1 graphene Landau level. 

PACS numbers: 73.43-f, 71.10.Pm 


I. INTRODUCTION 

The phenomenon of the fractional quantum Hall ef¬ 
fect (FQHE)P is the manifestation of one of the most 
remarkable many body states in nature. The electron’s 
spin degree of freedom adds to the richness of this phe¬ 
nomenon, and experiments have demonstrated the exis¬ 
tence of FQHE states with different spin polarizations as 
well as transitions between them as a function of the Zee- 
man energ 3 p®^. Analogous transitions have been f ound 
in multi-valley systems, such as AlAs quantum well d^^ * ^^ 
and two-dimensional electron system on an H-terminated 
Si(lll) surfac^^, where, in the SU(2) limit, the valley in¬ 
dex formally plays the same role as spin. 

The spin physics of EQHE is understood, qualitatively 
and semi-quantitatively, in terms of spinful composite 
fermions (CFsPt^. A composite fermiorPIJ^ is the 
bound state of an electron and an even number (2p) of 
quantized vortices. Composite fermions form Landau¬ 
like levels called A levels (ALs) in a reduced magnetic 
field given by 5* = 5 — 2pp(j)Q, where p is the electron 
density and = hcje is the flux quantum. The filling 
of composite fermions, z/*, is related to the electron fill¬ 
ing hy V = v"" jiflpy'' ± 1). The prominently observed 
EQHE at the the fractions v = sj(2ps ± 1) is explained 
as the z/* = s integer quantum Hall effect (IQHE) of 
composite fermions. Eor spinful composite fermions, we 
have s — where (s^) is the number of oc¬ 

cupied spin-up (spin-down) ALs. This immediately pre¬ 
dicts the allowed spin polarizations for any given fraction 
of this type, and the transition between two states with 


different spin polarizations occurs when CE ALs of oppo¬ 
site spins cross as the Zeeman energy is varied. Eor the 
V = s/(2s ± 1) EQHE states, modeling them as weakly 
interacting composite fermions correctly predicts the en¬ 
ergy ordering of the states, with the state with the small¬ 
est spin polarization having the lowest energy at small 
Zeeman energy and the fully spin polarized state win¬ 
ning at large Zeeman energies. 

There has been a resurgence of interest in EQHE 
in multicomponent systems as a result of the obser¬ 
vation of well developed EQHE in multivalley systems 
of graphen^^^H^. In particular, while the physics in 
the n = 0 graphene Landau Level (GLL) is consistent 
with the expectation from weakly interacting composite 
fermion^^, experiments by Amet et alW^ have indicated 
puzzling spin related behavior in the n = 1 GLL. These 
authors have seen extensive EQHE in the n = 1 GLL at 
the standard fractions of the form u = s/{2s ±1). In 
addition, they have detected a change in the transport 
behavior in tilted field experiments up to magnetic fields 
of B^45 T. Such tilted field experiments are interpreted 
in terms of a variation of the Zeeman energy, and the ex¬ 
perimental observations thus suggest that the spin degree 
of freedom remains relevant up to these magnetic fields. 
This is surprising for the following reason. The relevant 
parameter for the spin physics is /^ = Ez/(e^/e-^), where 
Ez is the Zeeman splitting, e is the di electric constant 
of the host material and I = ^JhcjeB is the magnetic 
length. In these experiments, the value of k at the high¬ 
est magnetic fields is hz ^ 0.07 (assuming a g factor of 
2.0 and e = 3.0). In GaAs, where the spin p hysics of 
the EQHE has been investigated in detaiP^®, the spin 
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degree has been found to be relevant only for k < O.Oi^. 
It is unclear why the spin degree should remain relevant 
up to much higher values of k in graphene than in GaAs. 
That was our primary motivation for undertaking the 
study reported in this article. 

Previous theoretical work has addressed this issue, but 
not conclusively. Using the methods of density matrix 
renormalization group (DMRG), Shibata and NomursP^ 
have studied the FQH states in n = 1 GLL in the torus 
and spherical geometries for systems up to 12 particles 
and surmised that the ground state at = 2/3 and 2/5 is 
fully polarized even at zero Zeeman energy. However, as 
we find below, the resolution to this issue is subtle and it 
requires much larger systems to conclusively determine 
the spin polarization of the ground state in the thermo¬ 
dynamic limit. Tdke and JairP^used the methods of GF 
theory in the spherical geometry, and concluded that the 
ground state in the n = 1 GLL is unpolarized at zero 
Zeeman energy and becomes fully polarized only at a fi¬ 
nite Zeeman energy. This work employed an interaction 
in the n = 0 LL that mimicked the Goulomb interaction 
in the n = 1 GLL (because the relevant wave functions 
are most easily evaluated within the n = 0 LL). However, 
in the course of our present study we have found that the 
interaction used in Ref. |34| was not sufficiently accurate 
for the purposes of spin physics, which is governed by 
very small energy differences. 

In this work we use a much more accurate interaction, 
and consider a wider range of filling factors than before. 
We also carry out a detailed comparison of the Moore- 
Read (MR) Pfafhan state and the GF Fermi sea (GFFS). 
We use methods similar to those in Refs. [34] and (35] 
Given that the observed fractions in the n = 1 GLL be¬ 
long to the standard u = sj (2ps±l) sequences, it appears 
very likely that the GF theory is qualitatively valid in the 
n = 1 GLL, but it is a priori unclear how accurate it is 
quantitatively, and to what extent the residual interac¬ 
tion between composite fermions may be neglected. 

Our conclusions, briefly, are as follows: 

We find that the u = s/{2s ± 1) FQHE states in the 
n = 1 GLL are fully polarized even at zero Zeeman en¬ 
ergy. The situation is somewhat tricky for both 2/3 and 
2/5, where the fully-spin-polarized and the spin-singlet 
states are almost degenerate insofar as their Goulomb 
interaction energy is concerned, and exact diagonaliza- 
tion studies are not able to decisively ascertain which 
of the two is the ground state. It is necessary to go to 
large systems using GF diagonalization (GFD) to deter¬ 
mine which of the two has lower Goulomb interaction 
energy. Our calculations thus strongly suggest that the 
tilted field dependence observed in Ref. [31] is unrelated 
to the spin physics of composite fermions (and as yet not 
understood). 

We further find that the fully spin polarized states 
are extremely accurately described in terms of composite 
fermions. The partially polarized states are also well de¬ 
scribed in terms of composite fermions. Taken together, 
these results imply that even though composite fermions 


are formed in the n = 1 GLL, the model of weakly 
interacting composite fermions does not remain valid. 
The exchange interaction between composite fermions is 
strong enough to fully spin polarize the states. Analogous 
physics was found previousl}^^ at fractions u = s/(4s + l) 
in the n = 0 LL. 

We also obtain the dispersion of the neutral excitons 
at = s/{2s±l). For the spin polarized excit on, th e dis¬ 
persion is very similar to that in the n = 0 LlJ ^^F^ I 
spin reversed exciton in the n = 1 GLL does not show 
any “sub-Zeeman energy” spin roton, as has been found 
for many FQHE states in the n = 0 LLP^®^. This is con¬ 
sistent with the above conclusion that the fully polarized 
state is stable against spin reversal. 

We have also considered filling factor 1 /2 in the n = 1 
GLL and found that the GEES remains stable to pairing. 
The MR Pfafhan state is not stabilized. This conclusion 
is consistent with previous worlP^®^. Eurt her more, we 
hnd that the Eermi sea also remains fully spin polarized 
even at zero Zeeman energy, in disagreement with the 
previous workP^^®!. This behavior is to be contrasted 
with the n = 0 LL where the GEES with minimal spin 
polarization has the lowest interaction energy. 

Given that the some of the competing states are almost 
degenerate, one may suspect if LL mixing can play a 
crucial role in determining the ground state. As discussed 
previously in Ref. |44l a reliable treatment of LL mixing 
is a complicated task, especially because the LL mixing 
parameter is fairly large (greater than 1) in graphene. 
Nonetheless, we have estimated the corrections and found 
that the EQHE states remain fully spin polarized even 
when LL mixing is incorporated. 

Throughout this work we shall neglect disorder, which 
may have less effect on these thermodynamic transitions 
than on single particle excitation gaps. Einite width ef¬ 
fect are known to be significant for the spin physics of 
EQH states in GaA^, but ought to be negligible in 
graphene. 

A discussion is in order regarding the validity of our 
use of the SU(2) model in the n = 1 GLL. Much work 
has been done to understand the lifting of various degen¬ 
eracies arising f rom t he spin and valley degree of free¬ 
dom in graphen^^^J^HEll^ principle, one can envisage 
the following three situations: 1) both the Zeeman and 
the valley-splitting terms are sufficiently small that the 
full SU(4) symmetry is preserved; 2) spin degeneracy is 
preserved while valley degeneracy is lifted or vice versa, 
whereby the appropriate symmetry would be SU(2); and 
3) spin and valley degeneracies are both lifted. Experi¬ 
ments carried out at high magnetic fields (values similar 
to the ones of Amet et find that the n = 0 GLL 

splits into four levels, thereby lifting both the spin and 
valley degeneracy, while in the n = 1 GLL, only the spin 
degeneracy is liftecP^. 

We assume in our work parameters such that the 
physics is well described within an SU(2) model, i.e., at 
most two of the four Landau bands are relevant. To the 
extent that the GE physics is valid, the neglect of SU(3) 
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and SU(4) symmetries is clearly valid for states in which 
composite fermions fill one or two A levels (e.g. 1/3, 2/5, 
2/3), because the state with s filled A levels cannot ex¬ 
ploit an S\J{k) symmetry for k > s. FQHE states at 
fillings factors with a numerator three or higher, such as 
3/5, 3/7, 4/7, 4/9, etc., can in principle support SU(3) 
and SU(^spin structureJ^. However, we have found in 
the pasfP^ that if a ferromagnetic state is preferred in 
a two-component system, then the inclusion of further 
components does not change this feature. Hence, given 
that we will find a ferromagnetic state even for the SU( 2 ) 
system in the n = 1 GLL, our neglect of SU(3) and SU(4) 
symmetries is a posteriori justified. (The only fraction 
in the n = 1 GLL where a non-ferromagnetic Goulomb 
ground state is feasible is at the filling fraction 2/7, but 
even here, a tiny Zeeman energy fully polarizes the spin, 
leaving only the SU( 2 ) valley degree to be considered 
here.) 

The principal motivation of our work are the exper¬ 
imental results of Amet et aM It is clear that the 
SU(3) and SU(4) physics are not relevant for these ex¬ 
periments for the following two reasons: (i) The exper¬ 
imentally observed FQHE states in the n = 1 GLL are 
precisely those expected from single component compos¬ 
ite fermions, namely s/{2ps ± 1). (ii) The evolution of 
transport gaps in the presence of a parallel magnetic field 
suggests that the spin degree of freedom plays a role even 
up to very high values of the magnetic field. It is therefore 
sufficient to consider SU( 2 ) physics for these experiments. 

We use the language of “spin” to label the two relevant 
indices, although our results below are also applicable if 
the two bands in question differ in their valley index (in 
that case “spins” would refer to the valleys, and “Zee- 
man splitting” to the splitting between valleys). At the 
end of the paper we briefly discuss the four component 
GEES state which involves the full SU(4) symmetry of the 
combined spin and valley degree of freedom of graphene. 
Recently states involving the full SU(3) or SU(4) sym¬ 
metry in graphene in the n = 0 GLL have been discussed 
in Ref. [51] 

We note that the physics of the n = 0 GLL, in the 
SU( 2 ) limit, is identical to that of n = 0 LL of non¬ 
relativist ic electrons in GaAs, to the extent that finite 
width and LL mixing effects may be neglected. Our con¬ 
clusions below for the n = 0 LL thus apply to both the 
n = 0 GLL and the lowest LL (LLL) of GaAs. 

The paper is organized as follows: In the next section 
we give a brief outline of the standard methods of exact 
and GF diagonlizations. We describe the effective inter¬ 
action used to simulate the n = 1 GLL in the n = 0 GLL. 
Sections [hi] and |IV| are dedicated, respectively, to the en¬ 
ergies of the ground state and low energy excitations at 
the filling factors u = s/{2ps^l). We include results from 
both GF theory and from exact diagonalization wherever 
these are available. In Section [V] we consider the ques¬ 
tion of stability of the GEES state against the paired MR 
Pfafiian state. We end with conclusions in Section ED 


II. EXACT AND CF DIAGONALIZATION 

All our calculations are carried out in the spherical 
geometr}!^ wherein N electrons reside on the surface of 
a sphere with a magnetic monopole of strength 2Qhc/e 
at its center. In this geometry the total orbital angular 
momentum L is a good quantum number and ground 
states are seen to be uniform and incompressible i.e., have 
L = 0 and have a finite gap to excitations. We consider 
states of all spin polarizations and apart from L we also 
use the total spin angular momentum quantum number 
S to characterize the states. The methods used in this 
work for exact and GF diagonalizati on are standard and 
details can be found in the literatur^^SHHH] briefly 
outline here some of the aspects that are relevant to the 
present work. 


A. Exact diagonalization 


In the spherical geometry, we consider N electrons con¬ 
fined to the surface of a sphere with the radial mag¬ 
netic field produced by a Dirac magnetic monopole of 
strength 2Q(j)o. The nth LL corresponds to the angu¬ 
lar momentum shell with / = Q + n (n = 0, 1, ...), 
which has 2 / + 1 states labeled by the z-component of 
the angular momentum m = —I, —/ + 1, ..., 1. The 
Goulomb interaction enters through the two-body ma¬ 
trix elements (mi,m 2 |E|ni 3 , 7714 ) connected to the pseu¬ 
dopotential V/ through the Clebsh-Gordan coefficients 
C'mima = {I, mi; I, m 2 \L, M = mi +m 2 ): 

(mi, m 2 V|m 3 , m 4 ) = ^ ( 1 ) 

In graphene the electron orbitals in all but the low¬ 
est LL are spinor states of the two-dimensional Dirac 
equation, with the two components corresponding to LL 
indices n — 1 and n: 


n, m)) 


V \n,m) ) 


( 2 ) 


The difference in LL indices on the right hand side com¬ 
plicates complicates the derivation of the {21 + l)-fold LL 
degeneracy from the model of a spherical surface pierced 
by the magnetic flux. Specifically, when Q defines a phys¬ 
ical monopole, then the two components have different I 
and hence different ranges of the allowed m’s. A common 
remedy has been to use pseudopot ential s V/ calculated 
(analytically) in the planar geometr}]^^^^^^!!] truncated to 
the allowed range /i <21. This will produce reliable re¬ 
sults for sufficiently large systems, but may have strong 
finite size corrections. We use an alternative approach 
following the direct analytical solution of the Dirac prob¬ 
lem on a sphere by JellaP^. Namely, we build the spinor 
wave function for graphene from a pair of orbitals with 
the same /, and thus with different Q. This leads to 
the expression for a spinor matrix element in graphene 
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through the scalar matrix elements in GaAs, all derived 
in spheri cal geo metry and hence all behaving properly at 
any rangd^^^ 

4 ((mi,m 2 || 4 ^||m 3 ,m 4 )) 

= (n — 0 , mi; n — 0 , m2|V^|n — 0 , m3; n — 0 , m4) 

+ (n — 1, mi; n — 0, m2|4^|n — 1, m3; n — 0, m^) 

+ (n — 0, mi; n — 1, m2|4^|n — 0, m3; n — 1, m^) 

+ (n — 1 , mi; n — 1 , m2|4^|n — 1 , m3; n — 1 , m^) . ( 3 ) 

We stress that all orbitals appearing in the above scalar 
matrix elements have the same /, which implies variable 
Q = / — nor / — n + 1, and that all scalar matrix elements 
are the spherical integrals. We first calculate the matrix 
element in units of e^/ei?, where R is the radius of the 
sphere. The matrix element is then expressed in the units 
^where the magnetic length t — Rj y/Q depends on 
the flux 2Q. We choose the mean of the two flux values to 
define the magnetic length = Rj^/Q^^ where is 

defined as the average of the two flux values, i.e., = 

(/ + (/ — l))/2 = / — 1/2 and express the energies in units 
of 

To take into account the presence of the positive back¬ 
ground charge density we simply subtract Ne^/2eR from 
the per-particle electron-electron energy (adequate for 
systems with zero widtlP^, where again we use R = 
Finally, to incorporate the density correction 
we rescale the energy by multiplying the backg round sub- 
tracted per-particle energy by the factor of yj2vQ^^/N. 
We note that this scheme produces correct pseudopoten¬ 
tials for large values of the relative angular momentum 
/i, an approximately linear dependence of the per-particle 
ground state energies on 1/A^, and identical thermody¬ 
namic energies as obtained using truncated planar pseu¬ 
dopotentials. 


B. Effective interaction 

The problem of interacting electrons in higher LLs is 
mathematically equivalent to a problem of electrons re¬ 
stricted to the LLL of GaAs (or the n = 0 GLL) but in¬ 
teracting via an effective interaction which repro¬ 

duces the higher LL Haldane pseudopotential^^. This 
turns out to be the most efficient way of doing calcula¬ 
tions within the GF theory, as the GF wave functions are 
most easily evaluated in the n = 0 GLL. In the spherical 
geometry the pseudopotentials depend on the size of the 
sphere. Gonsequently, one must find an effective interac¬ 
tion for each system size. To circumvent this problem we 
take the pseudopotentials from the disc geometry and use 
them for the spherical geometry. This procedure should 
give the correct results in the thermodynamic limit, and 
we carefully evaluate below the thermodynamic limits for 
all our GF calculations. The exact Coulomb pseudopo¬ 
tentials of graph ene o n the disc geometry are well known 
(see for exampl^^^^^^H^ . For the n = 1 LL these are 


Coefficient 

Value 

Bi 

1.0 

Bs 

0.5 

Bs 

0.5625 

Co 

42.55210599 

C 2 

-284.7958520 

C 4 

376.9514252 

a 

-179.3196140 

Cs 

36.58286249 

Cio 

-3.232277646 

Ci2 

0.1004240132 


TABLE I: Coefficients of Eq. which produces the effective 
interaction used to simulate the physics of the n = 1 GLL in 
the lowest LL. 


given by: 


T/(gr,l) 

* m 


rll (gr,l) 


Lr J 


m 



15m 153\ r(m-|) 

8 ^ 256y2r(m + l) ^ ^ 


The following form for the effective interaction in th e low- 
est LL is used to simulate the physics of n = 1 GLlJ^^^ 


W®(r) 


r 


Bs 

\/r^ 1 


^5 


k=6 

+E;C'2fer^''exp(-r2) 

k=0 


(5) 

where r is the distance in units of the magnetic length 
i. The coefficients H^’s and Q’s are listed in Table 
|Ij Hi, Bs and H 5 are fixed by the long range part of 
the Coulomb interaction while the C^’s are obtained 
by matching the first seven pseudopotentials of 
with those of the effective interaction. We find that 
fixing more than the first seven pseudopotentials leads 
to a highly oscillating form of the real space interaction 
which is undesirable, while fixing fewer pseudopotentials 
increases the error in the remaining pseudopotentials. 
Fig.E] shows the comparison of the pseudopotentials of 
this effective interaction in the lowest LL with those 
of the Coulomb pseudopotentials in the n = 1 GLL in 
the disc geometry. We see that the difference between 
the two is less than 0.06% for all values of m. We use 
this effective interaction to study large systems in the 
spherical geometry. This interaction is more accurate 
than the effective i nterac tion used in a previous study by 
some of the authorand the results below supersede 
the previous results for the n = 1 GLL. 


C. CF theory 

The CF theory states that the system of strongly in¬ 
teracting electrons in a magnetic field at a filling fac¬ 
tor u = s/{2ps ± 1 ) can be mapped into a system of 
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FIG. 1: (color online) Comparison between the exact pseu¬ 
dopotentials for the Coulomb interaction in the n = 1 GLL 
and the pseudopotentials of the effective interaction given in 
Eq.[5] in the n = 0 GLL in the disc geometry. 


weakly interacting composite fermions at a filling factor 
z/* = s = where (si) is the filling of spin-up 

(spin-down) composite fermions^. We denote this state 
by The ground state at these special filling fac¬ 

tors is described by the CF wave function: 

( 6 ) 

where J = ni<j</c<Ar(^i ~ ^k) is the so-called Jastrow 
factor, Zi denotes the coordinate of the electron, 
is the Slater determinant of s {s^ 
spin-up and spin-down) filled LLs for electrons 
and The states at z/ = s/{2ps — 1) 

require reverse flux attachment, and describe composite 
fermions in an effective field that points in a direction 
antiparallel to the applied magnetic field B. These 
wave functions involve LLL projection, for which three 
slightly different schemes have been employed. The most 
accurate one is the so called “hard core” projection, 
which has been shown to be very accurate^, but is 
not easy to implement. We will therefore not consider 
the hard-core projection here. Even for the non-hard¬ 
core projection, two m ethod s have been employed, 
namely direct proje ctior l^^ * ^^ * and the Jain-Kamilla 
projectHere, for technical reasons, we 
will exclusively consider the wave functions with Jain- 
Kamilla projection, which produces less accurate results. 
The details of the projectio n method have been outlined 
in the literatur j^^ l ^^ l ^^ l ^ ^^ and will not be repeated here. 

In addition to the above wave functions for the in¬ 
compressible states at z/ = s/{2ps ±1), we also study 
trial wave functions for the state at u = 1/2. For the 
CFFS, we will consider wave functions of the form in 
Eq.j^for which composite fermions experience zero fiuJ^ 
and then take the thermodynamic limit to obtain the 
energy of the state. We also study the Pfaffian wave 


functiorP^ Pf [M^j] J^, where Pf denotes the Pfaffian and 
Mij = {zi — Zj)~^^ which describes a chiral p-wave paired 
state of composite fermions,l^and is the most promising 
candidate for the EQHE state at u = 5/2 in GaAs sys¬ 
tems. 

Eor the GEES, we have considered for ff>(H* = 0) a 
single component state (which approaches in the ther¬ 
modynamic limit a fully spin polarized GEES), a two 
component state with equal number of particles in both 
components (spin singlet GEES) and a four component 
state with equal number of particles in all four compo¬ 
nents (SU(4) singlet GEES). In all cases, we work with 
filled shell states, so the state is uniform and isotropic. 

The total energies also include contributions from the 
positive neutralizing background, which we take into 
account under the assumption that a positive charge 
of Ne is uniformly distributed on the surface of the 
sphere. The net energy of N particles is given by: 
En = Eel-el — where the first term gives the 

contribution from electron-electron interaction, which we 
evaluate using the methods of GE theory, and the sec¬ 
ond term takes into account the electron-background and 
background-background interactions. To make a com¬ 
parison with experiments we need to find the energies of 
the ground states in the thermodynamic limit. The den¬ 
sity for a finite system depends on the number of par¬ 
ticles N and is different from its thermodynamic limit. 
Therefore, to reduce the N dependence we first apply 
a density correction to the energy, redefining it to be^ 
En = and then extrapolate to ^ oc. 

We note that the energies extrapolated to the thermody¬ 
namic limit with and without the density corrections are 
slightly different from each other. All numbers quoted in 
this work include density correction. 

Note that the background subtraction for exact diago- 
nalization is done slightly differently from the procedure 
mentioned here. Eor exact diagonalization the systems 
considered are not very large, so a more careful consid¬ 
eration is required while doing the background subtrac¬ 
tion. Eor results obtained using GE theory we go to much 
larger systems whereby the difference from these 0(1) 
corrections are negligible. Also since we take the ther¬ 
modynamic limit, the results are expected to be inde¬ 
pendent of these 0(1) corrections. We also note that the 
energies for a finite system obtained below from exact di¬ 
agonalization and GE theory are not directly comparable 
since the pseudopotentials of the effective interaction for 
a finite sphere are different from the n = 1 GLL Goulomb 
disk and spherical pseudopotentials used in the exact di¬ 
agonalization. The two results, however, should match 
in the thermodynamic limit. 

We evaluate energies of the wave functions given in 
Eq. using the Metropolis Monte Garlo method where 
in we sample with the above wave functions and use the 
yeff(^) interaction as the Hamiltonian. These energies 
are denoted by “GE w.f”. and are given in Tables [IH| and 


IV To improve on the above wave functions we use the 


method of GE diagonalization which is described next. 
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D. CF diagonalizatioiJ^ 

We take the V^^{r) interaction as the Hamiltonian and 
evaluate the energies of the states using the Metropolis 
Monte-Carlo method as follows: First, construct simul¬ 
taneous eigenstates of the and S‘^ operators in the 
corresponding IQHE system. These states are then mul¬ 
tiplied by the factor of and projected onto the low¬ 
est LL, a procedure known as composite-fermionization. 
This gives us the required (1/^,5'^) eigenstates since 
commutes with both these operators. The set of basis 
states are obtained by taking all possible (1/^,5'^) 

eigenstates. Evaluation of the Hamiltonian matrix in this 
basis involves multi-dimensional integrals for which we 
use the Metropolis Monte-Carlo method. In general, CE 
wave functions are not orthogonal to each other, therefore 
the Gram-Schmidt orthogonalization procedure is imple¬ 
mented to find the final eigenenergies. To improve the ef¬ 
ficiency of the computation we do the calculation within 
the subspace of (L^, S‘^) eigenstates. As was stated above 
the lowest LL projection is carried out using the Jain- 
Kamilla method details of which can be found in the 
literatur d^^ l ^^ l ^^ l 

The ground state wave functions can be improved by 
including CE excitons (where a single CE exciton is a 
pair of CE particle and hole) in the basis and perform¬ 
ing CEE^Sl. A single CE exciton does not couple to the 
ground state because its smallest angular momentum is 
one. Therefore, we need a minimum of two CE exci¬ 
tons to improve the ground state. In Eig. we show 
the excitations we considered to improve the 2/5 spin- 
singlet state. The Hilbert space grows very quickly with 
the number of excitons included in the basis for CED, 
so we restrict ourselves to at most two excitons. Among 
the ground state wave functions shown above, the fully 
polarized ones are extremely accurate, so this procedure 
of including two CE excitons in the basis of CED only 
marginally improves the ground state energy of the fully 
polarized state. Eor the partially polarized states the im¬ 
provement is more significant - for the spin-singlet states 
the energies improve by around 10%. (Note that the CE 
excitons in which the constituent CE particle and CE hole 
are separated by two ALs have L > 2 and therefore do 
not renormalize the incompressible ground state. Spin- 
flip excitons produce states in the degenerate {2S + 1) 
multiplet and therefore have identical ground state en¬ 
ergy but different Sz value. Consequently the two afor- 
mentioned excitons are not be included in the Hilbert 
space for CED.) We have performed CED for states at 
2/5 and 3/7. 


III. GROUND STATE SPIN POLARIZATION 

We have considered the spin physics at many frac¬ 
tions. Eor the sequences u = s/{2s ± 1), we have in¬ 
vestigated 2/5, 3/7, 4/9, 2/3, 3/5, 4/7; for the sequences 
u = s/(4s±l), we have studied 2/9, 3/13, 2/7, 3/11; and 
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FIG. 2: 2/5 spin-singlet states. Panel (a) shows the schematic 
of the the 2/5 spin-singlet CF state. At the first order approx¬ 
imation, composite fermion diagonalization mixes the state in 
(a) with the states shown in panels (b), (c) and (d) to obtain 
a new ground state with lower energy than the state in (a). 
Reproduced from Ref. [44l 


for the sequences u = s/(65 ± 1) we have studied 2/13, 
2/11, 3/17. As discussed below in more detail, all states 
at z/ = s/{2s ± 1) are fully spin polarized, although some 
non-fully spin polarized states are possible for the states 
of composite fermions carrying four or six vortices. 

A. Exact diagonalization 

We have performed exact diagonalization for up to 
N = 28 {N = 14) for fully polarized (spin-singlet) 2/3 
and N = 16 {N = 12) for fully polarized (spin-singlet) 
2/5. Similarly we have also carried out exact diagonaliza¬ 
tion up to A = 24 {N = 14) for fully polarized (partially 
polarized) 3/5 and N = 18 {N = 11) for fully polarized 
(partially polarized) 3/7. The results are shown in Eig.[^ 
and [H where we use both the truncated disk and the 
spherical pseudopotentials. Both pseudopotentials con¬ 
verge to the same thermodynamic limit, explicitly con¬ 
firming the validity of our procedure. We show here only 
the contribution of the Coulomb interaction energy; the 
Zeeman contribution can be incorporated separately. For 
the n = 0 GLL, the spin singlet state has lo wer Co ulomb 
energy, consistent with previous calcnlation d^'^ l ^^ l How¬ 
ever, in the n = 1 GLL, the fully spin polarized state 
and the spin singlet state are almost degenerate in their 
Coulomb interaction energy for both 2/3 and 2/5. This 
is a surprising result, but may appear somewhat more 
natural within the CF theory, where both of these states 
map into filling factor z/* = 2 of composite fermions. Be¬ 
cause there is no exact symmetry relating the energies of 
the spin singlet and fully spin polarized states, the near 
degeneracy in their interaction energy is accidental, and 
larger systems should select one of the two as the ther¬ 
modynamic ground state. Our calculations below based 
on the CF theory and also our results on the spin wave 
dispersion strongly point to the fully spin polarized state 
having slightly lower energy than the spin singlet state 
at both 2/3 and 2/5. The analysis of 3/5 and 3/7 (which 
both map to z/* = 3 of composite fermions) is similar. 
Exact diagonalization suggests that the fully polarized 
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FIG. 3: (color online) Thermodynamic extrapolation of the 
ground state energies, obtained by exact diagonalization, for 
the fully polarized and spin-singlet states in the n = 1 GLL 
at filling factors 2/3 (panel b) and 2/5 (panel d). Also shown 
for comparison are the corresponding results in the n = 0 
GLL in panels a (2/3) and c (2/5). We show results obtained 
from disk (red squares) and spherical pseudopotentials (blue 
circles). The background subtraction and density correction 
is carried out as described in Section EE Additionally, with 
the green crosses we show the expectation values of the n = 1 
GLL Goulomb interaction with respect to the exact n = 0 
GLL ground states; the latter are known to be nearly exactly 
described by the wave functions of GF theory given in Eq. 

All energies are quoted in units of e^/e^. 


and partially polarized states at these filling factors are 
almost degenerate, and the study of larger systems using 
the CF theory favors the fully polarized state. 

In Table [Ilj we show the overlaps between the ground 
states of the n = 0 and n = 1 GLL for different spin 
polarized states at various filling factors. The fully spin 
polarized states in the two Landau levels are nearly iden¬ 
tical. For the partially spin polarized or the spin singlet 
states, the overlaps between the n = 0 and n = 1 GLL 
states are not extremely high, but still high. We have 
also calculated the expectation value of the n = 1 GLL 
Coulomb interaction with respect to the exact ground 
state of the n = 0 GLL. These are shown as green crosses 
in Fig. [^and agree almost perfectly with the red squares. 
Because the n = 0 GLL states are nearly exactly cap¬ 
tured by the CF theory, these comparisons give us confi¬ 
dence in the applicability of the CF theory for the states 
in the n = 1 GLL as well. Note that this is very different 
from the behavior in GaAs, where the states in the n = 1 
LL are substantially different from those in the n = 0 LL. 


FIG. 4: (color online) Same as in Fig. |^but for z/ = 3/5 and 
ly = 3/7. 


B. CF wave functions 

The ground state energies obtained from CF wave 
functions of Eq. by linear extrapolation to the ther¬ 
modynamic limit of the different spin polarized states at 
various filling factors are shown in Tables |III| and |IV| 
Figs. and show both the linear and the quadratic ex¬ 
trapolations of the ground state energies. We note that 
the results obtained by the two procedures are consistent, 
except perhaps at u = 3/17 where the number of data 
points is too small for a definitive conclusion. In most 
cases, the energy difference is hardly affected by the ex¬ 
trapolation method, indicating that the slight curvature 
of the data points has no effect on the thermodynamic 
limit. We believe that this curvature is related to the 
background subtraction method, which does know about 
the effective n = 1 LL interaction. 

We find that for all states in the sequence u = s/{2s ± 
1), the fully polarized state has lower energy than the 
unpolarized state. In Fig. we plot the difference in 
the energies between the fully polarized and spin-singlet 
states at 2/5 as a function of 1/N. In the thermody¬ 
namic limit the fully polarized state is lower in energy 
than the spin-singlet state by ^ 0. 002e‘^/ei. As noted 
above, the wave functions of Eq. are very accurate for 
the fully spin polarized states but less so for non-fully 
spin polarized states, and one may ask if the energy or¬ 
dering may change in a more accurate calculation. In the 
next subsection, we present more accurate results from 
CF diagonalization, which show that the wave functions 
of Eq. overestimate the difference in energy between 
the fully polarized and partially polarized states but still 
produce the correct ordering. 

Next we discuss states of composite fermions carrying 
four and six vortices. Here again we find that the par¬ 
tially polarized states have higher energy as compared 
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ly 

Vt 

Ni 

2Q 

dimension 

K^s Vs^)l 

l(GVS^>l 



(2,0) 

2/5 

16 

0 

36 

155484150 

0.9998984626 

0.9999167952 

0.9998458746 

0.9998960866 

(1,1) 

2/5 

6 

6 

27 

2211680688 

0.9313810788 

0.9516235092^ 

U 

0.9965515722 

(2,0) 

2/3 

26 

0 

39 

259140928 

0.9998146190 

0.9998485003 

0.9998260141 

0.9998807588 

(1,1) 

2/3 

7 

7 

20 

280934870 

0.7882484953 

0.9522359035 

0.9999709236 

0.9279167113 

(3,0) 

3/7 

18 

0 

37 

386905330 

0.9999276937 

0.9999272643 

0.9997009362 

0.9997601209 

(2,1) 

3/7 

8 

3 

22 

17969272 

0.9133546180 

0.9374707945 

0.9998133677 

0.9981942320 

(3,0) 

3/5 

21 

0 

36 

155484150 

0.9998984626 

0.9999167952 

0.9998458746 

0.9998960866 

(2,1) 

3/5 

10 

4 

23 

383215178 

0.8480598539 

0.9209952153 

0.9998273715 

0.9819871145 


TABLE II: Table shows the Hilbert dimension and overlaps of the ground states of the n = 0 [superscript (0)] and n = 1 
[superscript (1)] GLL calculated using the spherical [subscript S] and truncated disc [subscript D] pseudopotentials for the 
largest systems considered in this work. ^For the 2/5 spin singlet state listed above, we do not have the vector We 

presume it is almost identical to and quote the numbers and 


to the fully polarized, with a sole exception at v = 2/7, 
where we find that the spin-singlet state has lower energy 
in comparison to the fully polarized state. We predict 
the the critical Zeeman energy for the transition from a 
spin-singlet state to a fully polarized state occurs at a 
Zeeman energy of 0.0014(4) ^|e^. We note that in the 
n = 0 GLL too, for composite fermions carrying more 
than two vortices, 2/7 is the oi^ filling factor in which 
the ground state is unpolarizecP^. As noted before, our 
calculation is inconclusive at v = 3/17; due to its mod¬ 
erate experimental relevance we do not pursue this issue 
further. 


C. CF diagonalization 


ly 

fully polarized (GF w. f.) 

spin singlet (GF w. f.) 

2/3 

-0.4122(4) 

-0.4097(7) 

2/5 

-0.3821(1) 

-0.3803(1) 

2/7 

-0.3551(1) 

-0.3558(1) 

2/9 

-0.32542(4) 

-0.32477(2) 

2/11 

-0.30176(1) 

-0.30170(3) 

2/13 

-0.281714(9) 

-0.28095(2) 


TABLE III: Thermodynamic limit of the n = 1 GLL energies 
obtained using the effective interaction of Eq. [^for the fully 
polarized and spin singlet states oXv — 2/(4p±l). All energies 
are in units of ^ |e^. These numbers are obtained using the 
GF wave functions in Eq.[^ Only for 2/7 does the spin singlet 
state have a lower energy than the fully polarized state. 


We have carried out the procedure of CFD to calculate 
the the ground state energies of the fully polarized and 
unpolarized states at two filling factors: 2/5 and 3/7. 
The linear extrapolation of the fully polarized and unpo¬ 
larized ground state energies obtained from CFD at 3/7 
are shown in Fig. The fully polarized and spin-singlet 
ground state energies are not well approximately by a 
linear function of 1/N for 2/5, but the difference in the 
two energies is. In Fig. ^we show this difference for 2/5. 
In the thermodynamic limit we find that the fully polar¬ 
ized state has lower energy than the spin-singlet state. 
For 3/7 the thermodynamic energy obtained from CFD 
for the fully polarized and partially polarized states are 
-0.3884(3) and -0.3877(1) respectively. As with 2/5, we 
again find that the fully polarized state has lower energy 
than the partially polarized state. The CFD results show, 
as anticipated, a smaller difference of ^ 0.001 ^ j be¬ 
tween the two energies than that obtained from the CF 
wave functions of Eq. 


D. Corrections from Landau level mixing 

The strength of LL mixing is quantified by a param¬ 
eter A, which is the ratio of the interaction to the ki- 


ly 

fully polarized (GF w. f.) 

partially polarized (GF w. f.) 

3/5 

-0.405(2) 

-0.399(2) 

3/7 

-0.3845(4) 

-0.3819(3) 

3/11 

-0.3489(3) 

-0.3488(2) 

3/13 

-0.32919(5) 

-0.32861(8) 

3/17 

-0.29803(9) 

-0.29751(7) 


TABLE IV: Same as in Table [Hi] but for z/ = 3/(6p ± 1). We 
find that the fully polarized state has lower energy than the 
partially polarized one for each hlling factor. 


netic energ}!^. In graphene this parameter is given by 
A = Iei)I{hvYIi) = e^/(fieuF), which also equals the 
graphene fine-structure constant. Unlike in the case of 
parabolic bands (as would be appropriate to semicon¬ 
ductor quantum wells such as GaAs), A is independent 
of the external magnetic field for electrons with a linear 
dispersion. For different substrates the values of A are 
as follows: A ~ 2.2 for suspended graphene, A ~ 0.9 for 
graphene on sili con d ioxide, A = 0.5 — 0.8 for graphene 
on boron nitridd^^^. 

We note that unlike in the n = 0 GLL, in n = 1 GLL, 
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FIG. 5: Extrapolation of the energy of the wave function in Eq. [^to the thermodynamic limit in the n = 1 Landau level of 
graphene for EQHE states in the sequence v — sl{2ps — 1). We show the results of both linear and quadratic hts. 



EIG. 6: Extrapolation of the energy of the wave function in Eq. [^to the thermodynamic limit in the n — 1 Landau level of 
graphene for states in the sequence v = s/{2ps + 1). We show the results of both linear and quadratic hts. 


LL mixing does produce an effective three body interac¬ 
tion in addition to the shifts to the pairwise interactiorP^. 
We have calculated the relevant (short range) pair and 
triplet amplitudes per electron in the considered ground 
states known for a series of system sizes then ex¬ 


trapolated these amplitudes to an inhnite system, and 
hnally convoluted the extrapolated amplitudes with pair 
and triplet pseudopotentials of the effective interaction 
which accounts for Landau level mixin^^. 

Fig. |9] shows the energy shifts induced by LL mixing 
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1/N 

FIG. 7: (color online) Extrapolation to the thermodynamic 
limit for the difference in the fully polarized and spin-singlet 
zy = 2/5 states in the n = 1 graphene LL. The red symbols 
show the energies obtained from the wave functions in Eq. 
while the blue symbols are obtained from CED. The energies 
are obtained using the effective interaction of Eq. The 
ground state is seen to be fully polarized. 



EIG. 8: (color online) Extrapolation to the thermodynamic 
limit for the fully polarized (blue) and partially polarized 
(red) V — 7 GED energies for the n — 1 GLL obtained 

using the effective interaction of Eq. The ground state is 
seen to be fully polarized. 


for the fully polarized and unpolarized ground states in 
n = 1 GLL. As is evident from the figure, LL mixing 
favors fully polarized states over unpolarized states. We 
note that Ref. iTOl does not give the 3-body pseudopoten¬ 
tial for spin s = 1/2 and m = 3, which we arbitrarily set 
to zero. However, the corresponding amplitudes are very 
small, so we do not expect that adding the actual value of 
V{s = 1/2, m = 3) would affect the result. We end this 
subsection with the caveat that the corrections in Ref. |70| 
are obtained within a perturbative scheme, and we have 
not tested to what extent this approach remains applica¬ 
ble for the realistic values of the LL mixing parameter A 



EIG. 9: (color online) Thermodynamic extrapolation of the 
exact Goulomb ground state energies per unit LL mixing pa¬ 
rameter A for the fully polarized (cirlces) and unpolarized 
states (squares) in the n — 1 GLL at filling factors 2/3 and 
2/5 (a) and 3/5 and 3/7 (b). 


(See Ref. |44|for a detailed discussion). 


IV. EXCITATIONS 

The FQHE ground state at z/ = s/(2ps ± 1) are ex¬ 
plained as s ALs of composite fermions carrying 2p vor¬ 
tices. Furthermore, their low-energy excitations are ob- 
tained by exciting composite fermions across AL d^^ * ^^ ^. 
Our motivation for investigating the excitations comes 
from the fact that spin reversed excitations of composite 
fermions can reveal an instability of the fully polarized 
ground state at sufficiently small Zeeman energies. Such 
an instability is signaled by the presence of a sub-Zeeman 
energy rot on in the dispersion of the spin reversed neu¬ 
tral excitoiP^® theoretically, it is necessary to dress the 
single spin reversed CF exciton by other excitations to 
capture the physics of the spin roton. 

We have studied two kinds of low-energy collective 
modes at the filling factors u = s/{2s ± 1) assuming the 
ground state to be fully polarized. These are: 

(i) Spin-conserving: a CF is excited from the {s — l)t AL 
to the st AL. These excitons form the magnetoro- 
ton/magnetoplasmon mode. Their energies have been 
calculated theoretically for many of the fractions in the 
LLlP. Resonant inelastic light scattering experiments 
have observed these modes and their measur ed energ ies 
are consistent with the theoretical predict iond^^MIll, 

(ii) Spin-reversing: a CF is excited from the {s — l)t AL 
to the 0| AL. These excitons form the spin-flip mode. 
The energy of these modes are evaluated numerically 
in the lowest LL and it was found that for n > 2 this 
excitation has a lower energy than the ground state 
up to a finite Zeeman energy showing that the fully 
polarized ground state is unstable to these excitations^. 
This is consistent with the fact the ground state i n the 
LLL at zero Zeeman energies is spin unpolarizecP^^^. 
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FIG. 10: (color online) Spin wave (black hexagram), spin- 
conserving (blue squares) and spin-reversing (green penta¬ 
gram) collective mode energies for A/’ = 30atz/ = 2/5 for 
the n = 1 GLL. These energies are obtained using the wave 
function in Eq. and the effective interaction of Eq. in (top 
panel). Eor comparison the corresponding modes in the lowest 
LL are shown in the bottom panel. 


For comparison we have also evaluated the energy of the 
spin-wave mode. 

In Fig. [^we show the collective modes in the n = 1 
GLL for a fully polarized at = 2/5 obtained using the 
CF theory with the effective interaction of Eq. In 
the bottom panel of the same figure, the corresponding 
modes in the lowest LL are shown for comparison. In 
contrast to the lowest LiP, the spin-flip modes in the 
n = 1 GLL does not have a roton minimum below its long 
wave length energy. Consequently, unlike the lowest LL, 
there is no spin-flip instability in the second GLL. It is 
interesting to note that in the n = 1 GLL the large wave 
vector gaps for the spin conserving and spin reversing 
modes are almost identical, in contrast to the LLL where 
the two gaps are significantly different. This is seen in 
the spectrum shown in Fig. p!Q| where the two collective 
modes are seen to approach each other in the n = 1 GLL 
in the large wave vector limit. 

We have also obtained the spin wave dispersion for 



FIG. 11: (color online) Gomparison of the spin wave disper¬ 
sions in the LLL (left panels) and n = 1 GLL (right panels) at 
various filling factors in the sequence u = s/(2s± 1), obtained 
from exact diagonalization. The spin-wave in the LLL shows 
a spin-flip instability for s > 2, while the n = 1 GLL supports 
a robust spin wave. 


the fully polarized state in the n = 0 and n = 1 GLL 
from exact diagonalization, shown in Figs. and 
It is evident that the fully spin polarized state at filling 
factors u = s/{2s ± 1) for s > 2 in the n = 0 GLL is 
unstable to a spin-flip instability at vanishing Zeeman 
energies whereas that in the n = 1 GLL remains stable. 


V. IS THERE A PAIRING INSTABILITY OF 

COMPOSITE FERMIONS IN THE n = 1 GLL? 

Given that filling factor ly = 1/2 in the n = 1 LL 
of GaAs produces an incompressible state, it is natu¬ 
ral to ask if the same is true in graphene. This ques¬ 
tion has been addressed in the past using both exact 
diagonalizatioiP^ and methods of CF theor}®, and it 
was concluded that the Pfaffian states is not favored. For 
completeness, we have addressed this issue with our more 
accurate real-space interaction. 

Figj^ shows an extrapolation of the energies of the 
CFFS states (with different spin / pseudospin polariza¬ 
tions) and the MR Pfaffian state. The Pfaffian wave 
function is clearly ruled out. Furthermore, the fully spin 
polarized CFFS has the lowest energy. Even though the 
CF wave function for the spin singlet CFFS is not as ac¬ 
curate as that for the fully polarized CFFS, the energy 
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FIG. 12: (color online) Comparison of the spin wave disper¬ 
sions in the LLL (hlled symbols) and n = 1 GLL (open sym¬ 
bols) at various filling factors in the sequence v — sji^s zb 1) 
for different system sizes. The dispersions are obtained from 
exact diagonalization. 


difference between them is sufficiently large that we are 
confident in concluding that the ground state at z/ = 1 /2 
in the n = 0 GLL is t he ful ly spin polarized CFFS even 
at zero Zeeman energ}^^^. For contrast, we also show 
results for two different interactions. One corresponds 
to the n = 0 LL, where also the Pfafhan has very high 
energy, but the ordering of the CFFS states is reversed 
(with the least polarized state having the lowest energy). 
The second interaction is that of n = 1 LL of GaAs, 
where the Pfafhan wave function is seen to produce the 
lowest energy among the trial states studied, consistent 
with a previous stud^P^. 

In our studies of z/ = 1/2 we have neglected the effect 
of LL mixing. Ref. |70] included the effects of LL mixing 
to the lowes t ord er and found that it favors the anti- 
Pfaffian stat eJ^^ * ^^ (particle-hole conjugate of the Pfaffian 
state) over the Pfaffian state in the n = 1 GLL for a 
certain range of the LL mixing. The question of whether 
LL mixing can drive the paired state below the fully spin 
polarized CF FS ha s not been explored. 

TheoreticaP^^^ and experimentaP^ evidence exists for 
the formation of the Pfaffian state in hilayer graphene. 
We have focused only on monolayer graphene here. 


VI. CONCLUSIONS 

We have used methods of exact diagonalization and CF 
theory to obtained the ground state energies of differently 
spin polarized FQHE ground states along the sequence 
z/ = s/{2s ± 1) in n = 1 GLL, as well as their neu¬ 


tral excitations. Our principal conclusion is that these 
states are fully polarized even at vanishing Zeeman en¬ 
ergy. This is in stark contrast to the n = 0 LL, where, at 
very small Zeeman energies, the CF state with the small¬ 
est polarization is the ground state. We further establish 
that these states are well described by the CF theory, 
but the model of weakly interacting composite fermions 
breaks down in n = 1 GLL. In the absence of their in¬ 
teraction, composite fermions would occupy the lowest 
ALs. When multiple ALs are filled, the ground state 
in the absence of Zeeman splitting will be the one with 
the lowest spin polarization. Hence, the calculated full 
spin polarization of the actual ground state of electrons 
with Coulomb interactions implies sufficiently strong in¬ 
teraction among the composite fermions to overcome the 
single-CF AL splitting (i.e., the effective CF cyclotron 
energy). Furthermore, in the case of filled ALs, the rel¬ 
evant interaction among the composite fermions is their 
exchange. Therefore our calculations indicate that the 
exchange interaction between composite fermions is suf¬ 
ficiently strong to drive the system ferromagnetic even in 
the absence of Zeeman coupling. We have also estimated 
the corrections arising from LL mixing. We find that our 
conclusions regarding the spin polarization of the FQHE 
ground state remains unchanged even after incorporating 
the effect of LL mixing. 

We have considered two neutral excitations of the fully 
polarized ground state at u = s/{2s ±1): namely the 
spin-conserving exciton and the spin-flip exciton. We find 
that the spin-conserving exciton has a dispersion similar 
to that in the n = 0 LL. On the other hand, the spin-flip 
exciton in the n = 1 GLL does not show any spin-roton 
with a negative Coulomb energy unlike its counterpart 
in the n = 0 LL. This further corroborates the tendency 
of the ground state to be fully spin polarized at z^ = 
s/{2s±l). 

Our study further deepens the puzzle surrounding the 
experimental results of Amet et. ai^ discussed in the 
Introduction section. They found that the experimental 
data are sensitive to the Zeeman energy even at very 
large Zeeman energies, which led them to conclude that 
the spin was playing a role even at very high magnetic 
fields. Our study, on the other hand, indicates that the 
system is fully polarized, at least at the special fractions, 
even at zero Zeeman energy. An explanation of the Amet 
et al result remains elusive. 

Einally we looked at the ground state at half filling in 
the n = 1 GLL. Here, we find that CEES remains stable, 
and CE pairing is not favored. Eurthermore, the CEES 
remains fully spin polarized even at zero Zeeman energy. 
These facts should be contrasted with the n = 1 LL of 
non-relativistic electrons of conventional semiconductors 
where the CEES is unstable to a Pffafian type pairing, 
and the n = 0 LL of conventional semiconductors where 
the CEES is unpolarized in the absence of Zeeman cou¬ 
pling. 
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FIG. 13: Extrapolation of the ground state energy to the thermodynamic limit of the states at half filling (a) in the lowest 
Landau level, (b) in the second (n = 1) Landau level of GaAs, (c) in the n = 1 Landau level of graphene. The wave function in 
Eq. [^has been used to obtain the energies. The Zeeman energy has been omitted, and the density correction has been applied. 
The solid lines correspond to a linear fit, and the dashed lines to a quadratic one. In (b) the energies in the thermodynamic 
limit may differ from those reported previously in the literature because of different methods for background subtraction, but 
the energy differences are captured properly. 
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